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Abstract. In this paper we solve one dimensional trapped SU(2) bosons with 
repulsive ^-function interaction by means of Bethe-ansatz method. The features 
of ground state and low-lying excited states are studied by numerical and analytic 
methods. We show that the ground state is an isospin "ferromagnetic" state 
which differs from spin-1/2 fermions system. There exist three quasi-particles in 
the excitation spectra, and both holon-antiholon and holon-isospinon excitations 
are gapless for large systems. The thermodynamics equilibrium of the system at 
finite temperature is studied by thermodynamic Bethe ansatz. The thermodynamic 
quantities, such as specific heat etc. are obtained for the case of strong coupling limit. 
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1. Introduction 

In the past few years, Bose Einstein condensation (BEC) has witnessed a sequence of 
exhilarating experimental achievements, and much attention has been paid to the study 
of Bose systems. The possibility of BEC in one-dimension has been discussed for the 
non- interacting Bose gas [||, The role of dimensionality has been examined for the 
ideal Bose gas It is known the interaction between bosons plays an essential role 
in one dimension due to the strong constraint in phase space |Q. The Luttinger liquid 
properties of trapped, interacting quasi-one-dimensional Bose gas was discussed 0. 
It was shown under suitable experimental conditions that the system can be described 
as a Luttinger liquid and the correlation function of the bosons decays algebraically 
which prevents Bose-Einstein condensation. Recently, a two component Bose gas was 
produced in magnetically trapped ^^Rb by rotating the two hyperfine states into each 
other with the help of slightly detuned Rabi oscillation field 0, §]. It was noticed 
that the properties of the Bose system can be different from the traditional scalar Bose 
system once it acquires internal degree of freedom. The two-component Bose system 
on a circle (periodic boundary conditions) has been studied by means of Bethe-ansatz 
method ||T^. It was shown that the ground state is an isospin-ferromagnetic state which 
differs from the one-dimensional spin-1/2 Fermi system whose ground state is SU(2) 
singlet. 

In this paper, we solve trapped two-component bosons with 5-function interaction 
in one dimension by Bethe-ansatz method. On the basis of Bethe-ansatz equation, we 
discuss the ground state, low-lying excited states and the thermodynamics of the system 
at finite temperature, where some thermal quantities at low temperature are obtained 
explicitly. Our paper is organized as follows: In the next section we introduce the model 
and derive a set of non-linear equations for charge rapidity and isospin rapidity. In Sec. 
^, we explicitly show that the ground state is an isospin "ferromagnetic" state and 
manifest how the quantum numbers in Bethe-ansatz equation should be taken for the 
ground state. In Sec. We study the low-lying excited states extensively by analyzing 
the possible variations in the sequence of quantum numbers. Numerical results of energy- 
"momentum" spectra for each excitation are given. In Sec. |^ we discuss the general 



thermodynamics of the system by means of thermodynamic Bethe ansatz [|1^]. In Sec. 
^ we derive the free energy and specific heat in the case of strong coupling limit. In the 
last section a brief summary is given. 

2. The model and its Bethe-ansatz solution 

The Hamiltonian for a two-component bosons trapped in a potential well of infinite 
depth reads 

N r.2 N N 

1=1 * i=l i>j=l 
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where 

The Hamiltonian is not invariant under translation due to the presence of trapping 
potential, and the total momentum of the system is not conserved. However, the system 
is still invariant under the action of the permutation group Sn which makes it possible 
to employ the coordinate Bethe-ansatz approach. 

In the domain with Xi ^ Xj{i ^ j) and inside the potential well, the Hamiltonian 
reduces to the one for free bosons and it's eigenfunctions are therefore just the 
superpositions of plane waves. When two particles collide with each other at the same 
point, a scattering process happens. The Bethe-ansatz embodies that this process 
is purely elastic and the occurrence during the process is merely that the particles 
exchange their momenta. So for a given momentum k = [ki, k2, ■ ■ ■ , k^), the scattering 
momenta include all permutations of the components of k. Moreover, since the total 
momentum is not conserved in the present model, we will have more nondiffr active 
scattering momenta, e.g., for = 2, we have (/ci, ^2) {—ki, ^2), (^1, —^2), (— ^i, —^2), 
(/c2, ki), (— ^2, ki), (^2, —ki), or (— /c2, — ^i)- All the eight states correspond to the same 
energy kf + Hence, besides all possible permutations of the components of k, the 
scattering must include all possibilities of sign changes in the components of k. 

Now we consider the case of bosons. Considering the fact that the total 



momentum is not conserved, we adopt the following Bethe-ansatz form [14 



M^) = J2 MP, Q)e*(^'i^"^ X e CiQ) (3) 

where a = (ai,a2, . . . ,aAr), aj denotes the isospin label of the jth particle; Pk stands 
for the image of a given k := {ki, ^2, ■ ■ ■ , k^) by a mapping P G Wb and the coefficients 
A{P, Q) are functionals on Wb ® Wa- Here, Wb and Wa stand for the Weyl group 
of Bj\f and Aj^j^i Lie algebras respectively. The later is isomorphic to the permutation 
group Sn and the former consists of Sn and all possible sign changes. We emphasize 
that the sum runs over the Weyl group of Lie algebra Bn but the wave function is 
defined on various Weyl chambers C{Q) corresponding to the Weyl group of An-i Lie 
algebra. This is different from the situation of periodic boundary conditions. 

For a Bose system, the wave function is supposed to be symmetric under any 
permutation of both coordinates and isospin indices, i.e. 

{a'^)a{x) = M^) (4) 

Here {o'^4')a is well defined by ip^ja{o'^x), therefore both sides of (0) can be written out 
by using d). Furthermore, using the evident identity (Pkla'-Qx) = {a''Pk\Qx) and the 
rearrangement theorem of group theory, we obtain the following consequence from (^) 

A,{P,aV)=A^,,{a^P,Q). (5) 
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The 5-function term in the Hamiltonian (|1]) contributes a boundary conditions at 
hyperplane [a is a root of Lie algebra A^r-i), namely a discontinuity of the derivatives 
of the wave function along the normal of a Weyl hyperplane: 

lim [a ■ ViJa{x(a) + ea) - a ■ V^/'a(x(a) - ea)] = 2c?/'a(a;(„)) (6) 

where X(^a) ^ Pa and V := J2iLi ^{{d/dxi). Substituting @ into @, we find that 

z[(Pfc), - (PA;),+i] [Aa{P, a,Q) - A,{a,P, a,Q) - A,{P, Q) + A,(a,P, Q)] 

= c[Aa{P, Q) + Aa{a,P, Q) + Aa{P, Q) + Aa{a,P, Q)] (7) 

By making use of the relation and the continuity relation of wave function across 
{Qx)i = {Qx)i+i, we can obtain the following relations 

A,{a,P, Q) = SlAPk)AAP, Q) (8) 

^'^"^'^ ^" C-.[(Pfc).-(PfcM 

where Pa,a' stands for the matrix elements of the spinor representation of the 
permutation group. The relation provides for the coefficients A a relation between 
different Weyl chambers. Eq.(^) provides a connection between those coefficients which 
are related via any element of Weyl group Wb in the same Weyl chamber. 

The basic elements of the Weyl group Wb obey aV* = 1 and o"V*^^cr* = 
cr*+^crV*+^ as identities. These identities imply that Aa{(T^(T'^P,Q) = Aa{P,Q) and 
Aa(crV*+V*P, Q) = Aa((T*+VV*+^P, Q). Using (H) repeatedly, one can obtain the 
following relations: 

S'{(T'Pk)S'{ak) = I 

S\a'+^a'Pk)S'+^{a'Pk)S\Pk) = S'+\a'a'+^Pk)S%a'+^Pk)S'+\Pk) 

(10) 

where we have adopted the conventions S = mat{Sab), S"* = S" ® /, 5"*"*"^ = I ^ S (/is 
a 2 X 2 unit matrix). These relations are consistency conditions for the S-matrix. The 
second relation is called the Yang-Baxter equation. The concrete S-matrix in (^ fulfills 
these relations. 

Due to the infinite depth of the potential well, the wave function must vanish at 
the ends of the well, 

MiQ^)i = -L/2) = 0, MiQ^)N = L/2) = (11) 
These boundary conditions give rise to 

Aair'P, Q) = -e-^(^'=)iM„(P, Q) (12) 
AaiiP, Q) = -e^(^'=)-M„(P, Q) (13) 

where The relation ( |12|) together with eq.(|D determine 

the amplitudes A's up to an overall factor since a^, . . . , o"^^^ and generate the 
whole Wb- There are two further consistency conditions. One comes from the identity 
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{r^a^T^a^ = a^T^a^r^) of the Weyl group Wb leading to a reflection related Yang-Baxter 
equation which is fulfilled automatically due to the fact that the reflection matrix is just 
the unit matrix multiplied by a scalar function. The other one comes from eq.(^) which 
leads to an eigenvalue equation for the products of the S-matrices. This eigenequation 
can be diagonalized by means of standard quantum inverse scattering method [p!5[]. One 
finally obtains the following Bethe-ansatz equations 

i2kjL _ TJ ~ + kj + ki + ic kj - Xfi- ic/2 kj + - ic/2 
kj — ki — ic kj + ki — ic j^J-^ kj — + ic/2 kj + A^ + ic/2 

^ _ yr - ki - ic/2 + ki - ic/2 yr X^ - X^ + ic X^ + X^ + ic 

j-J^ X-y — ki-\- ic/2 X^ + ki + ic/2 A^ — A^ — ic A-y + A^ — ic 

where M denotes the total number of down isospins and A denote isospin rapidities 
which arise from the diagonalization conditions of quantum inverse scattering method. 
Taking the logarithm of eqs. (|14D gives rise to the secular equations: 



7r 



tan-^ + tan- f^^l - V ftan-^ + tan^ + 



1=1 -I -< M7^7 



c/2 / V c/2 



(15) 

where Ij is the quantum number for the charge rapidity kj and for the isospin 
rapidity. Concerning the property of the logarithm function, Ij and take integer 
values regardless of either the total number of particles or that of isospin down. Once 
the roots are solved from eqs. (|T3|) for a given set of quantum numbers {Ij, J^}, the 
energy can be calculated by = Xl^i ^f- 

3. The ground state 

For any set of quantum numbers {/j, J^} with the solution {/cj,A^}, the replacement 
of either Ij —Ij, kj —kj or —Jj, X^ —X^ keeps the secular equation 

(p!5|) invariant. So it makes no change to the energy, then we only need to consider 
the case of positive integer [|1^ for Is and Js. In the weak coupling limit c — 0, 
tan~-'^(x/c) 7rsgn(x)/2, eqs. (|T5D become 

N M 

nlj = kjL + - ^[sgn{kj - k/) + sgn{kj + ki)] - - ^[sgn(A;j - A^) + sgn(A;j + A^)] 

N M 

2^7 = X][^Sn(A^ - h) + sgn(A^ + ki)] - ^[sgn(A^ - A^) + sgn(A^ + A^)] (16) 

1=1 
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We may choose both the subscript of the Jj and in an increasing order so that the 
second equation of eqs. (|T6|) becomes 

N 

J2 sgn(A^ -ki) = 2J^ + 2-f - N -2 (17) 
1=1 

which also gives rise to 

N 

J][sgn(A-,+i - ki) - sgn(A^ - k)] = 2(^1 -J^) + 2 (18) 
1=1 

Thus, for J-y+i — = m, there must exist exactly m + 1 momenta of ki satisfying 
A^ <i ki <i A^-|_i. 

Similarly, from the first equation of eqs. (|T^) we have 

2L ^ 
2(^i+i-^i) {kj+i-kj)-2 = - J^[sgn(/cj+i-A^)-sgn(A:j-A^)](19) 

/i=l 

For Jj_|_i — Ij = n, there will be kj^i — kj = {n + m — l)'/r/L if there are m A^ satisfying 
kj < A^ < kj^i. Therefore, an isospin rapidity of value A^ always repels the quasi- 
momentum away from that value. The ground state for c — > does not allow any possible 
existence of A^ between ki and kjy. Thus from eq. ([T7| ) and eq. (|^) we can conclude 
that the total number of isospin down should be zero. The quantum numbers {Ij} of 
the ground state are given by the successive positive numbers /° = 1, Ij_^_i — Ij = 1. 
In the thermodynamic limit, as the difference between two adjacent fc's is small, the 
contribution of any existing A to the density of /c-rapidity at points k = X will brings 
about a rift (see fig. |1] right). Therefore, an existing A^ will suppress the density of state 
in fc-space at the point k = The more isospin rapidities there are, the higher the 
energy will be. In contrast to the spin-1/2 Fermi system whose ground state is SU(2) 
singlet, the ground state of trapped 2-component bosons is an isospin "ferromagnetic" 
state which is described by the quantum number, 

{J,}:={l,---,iV}, {J^} = empty (20) 

Since the existence of open boundary conditions makes the backwards scattering k ~k 
possible, there exists an additional term tan~^(2A;/c) which contributes a negative value 
to p{k), especially p(0) is obviously suppressed. This phenomenon is different from that 
for the system with periodic conditions. We interpret it being due to the confinement 
by the infinite-depth well, for example, the ground state for a single particle is not = 
but = tt/L in open boundary conditions. The density of states of the ground state for 
various coupling is plotted in fig.|l| the left. 
By introducing 



kj + kj+l^ 


1 






kj^i - 


kj 


' A^ + A^_(_i ^ 






2 y 




A-y 



Lp 

Lcr{'-^^^-^)=-—— (21) 
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the density corresponding to the quantum numbers (|2y) satisfies the integral equation 
Poik) = -- yKy^ik) + K,{k\k')po{k')dk' (22) 

71" ^ Jo 



where 



and 



KJx) = -- 



Knix\y) = Kn{x -y) + Kn{x + y) 



PoiQo are respectively the density and integration limit for the ground state. The 
Ki/2{k) term in ( P^ comes from the fact that tan~-^( ^^^^' ) but tan"-^'' ^^"'"^' '' ^ 

limit. The particle number per length is determined by 



c 



tan (^) when taking account of the term of / = j in the process of thermodynamic 



N 

D = -= I poik)dk (23) 







eqs. ( p2D and eq. ([23|) determine Qo and po{k). Then the energy can be calculated by 
^" epik)dk (24) 



L 







1 4 

which is explicitly -n'^D^il D) in the strong coupling limit c ^ 1. In the general 

3 c 

case one needs to solve the equations numerically. We show the ground state energy for 
various densities D = 1.0, 0.75, 0.5 in fig.^ 



4. Low-lying excited states 

The excited states are obtained by the variation of the configuration {Ij, J^} from that 
of the ground state. The simplest case is to remove one of Is from the configuration of 
the ground state and add a new outside the original sequence. We call this holon- 
antiholon excitation which is described by 

{/,} = {l,---,n-l,n + l,...,Ar,/4 

with |/„| > N and M = 0. To investigate the excited states, we shall consider systems 
of finite size first, then take thermodynamics limits. Although the total momentum 
kj is no more a constant in open boundary conditions, Ylj \kj\ is still a constant for 
a given energy eigenstate. We plotted the numerical results for energy- "momentum" 
spectrum in fig.^ and ^ for = 40, L = 40 where the x-axis represents the quantitive 
change ^j(|fcj| — |^j|)- We notice that the overall structure of the spectrum does not 
change substantially between weak (c = 1) and strong (c = 10) coupling regimes. For a 
finite-size system, the gap of holon-antiholon excitation opens, and it's dependence on 
the inter-particle interaction is shown by the upper line of fig.^. 

In the thermodynamic limit, it is plausible to calculate the excitation energy by 
making p{k) = po{k) + pi{k)/L where poik) is the density of the ground state. By 
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creating a hole inside the quasi Fermi sea k G [0, Qo] and an additional kp > Qq outside 
it, we have 

pi{k) + 5ik-k) = dk' pi{k')Ki{k\k') + Ki{k\kp) (25) 







The excitation energy consists of two terms /\E = f k^p{k)dk — f k'^po{k)dk = 
£h(k) + Saikp). The holon energy Eh and antiholon energy Saikp) = —Shikp) are given 

by 



Jo 

ef,{k) = -P+ / k^p'l{k,k)dk, 
Jo 

p'lik, k) = -Ki{k\k) + / ° Ki{k\k')p\{k'\k) (26) 

Jo 

Another interesting excitation is to flip one isospin down, i.e., M = 1. As discussed 
above, a single A can suppress the density of state at = A. We plot the density of 
the lowest energy state for various couplings in flg.|I| the right with = 100, L = 100 
and Ji = 30. The holon-isospinon excitation is characterized by the Young tableau 
[A^ — 1, 1], accordingly, J's take A^ distinct integers and Ii < Ji < In, namely 

= -iV/2 + (l<ji<Ar + l) 

I, = Ij.i + 1 + 6j^j, {j=2,...,N) 

where 6a,i3 is the usual Kronecker symbol. The excitation spectra are plotted in flg.§, 
for L = 20, A^ = 20, c = 10 and c = 1 respectively. The lowest energy state of one 
isospin down is to make the absolute value of Ji as large as possible in order to avoid 
the enhancement of energy caused by the suppression of the fc-distribution. In the 
thermodynamic limit, it can be shown that this mode is gap less. For a flnite particle 
system, however, it still has a gap of order 1/L which decreases more quickly than that 
for periodic conditions for particle number ranging from small to large. Figure |^ shows 
it's dependence of the gap on the interaction. We presented in flg.^ the behavior of the 
finite-size spin gap as a function of 1/L at c = 1.0, 10.0, 100.0. 

In the thermodynamic limit, we should take into account of the hole in the /c-sector 
and the pi{k) should be solved from 

pi{k) +5{k-k)= [ Ki{k\k')pi{k')dk' - Ky2ik\X) (27) 



One can find that the energy of holon-isospinon excitation consists of two terms 
AE = Shik) + eiso{X). The Sh is determined by eqs. ( ^61) and the Siso by Siso^X) = 
J k^p\'°{k, \)dk with 



Jo 



pr{k,X) = -Ky2{k\X)+ K,{k\k')pr{k\X)dk' (2^ 
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5. Thermodynamics at finite temperature 



For the ground state (i.e. at zero temperature), the charge rapidity fc's are real roots 
of the Bethe-ansatz equations ([T^). For the excited state, however, the isospin rapidity 
can be complex roots |Tl| which always form a "bound state" with several other A's. 
This arises from the consistency of both sides of the Bethe-ansatz equations |jl2| . The 
n-string of A rapidity is defined as 

= a: + (n + 1 - 2j)ic/2 + 0{exp{~6N)) (29) 

where j = 1,2 ■ ■ -n. The number of total down isospin is determined by 



(30) 



n=l 



where M„ denotes the number of n-strings. The eqs. (|T^ become 



N 



tan 



-1 1 



+ tan 



-E 

an 

ttJ? = tan"^ 



tan 
A" 



nc/2 



+ tan 



nc/2 



nc/2 



N 



mbt 



mnt 



tan 



1=1 



tan 
A" - A 



1 / ^'a 



tc/2 



nc/2 



+ tan^ 



tan 

1 I + \ 

tc/2 



K + k 



nc/2 



(31) 



where 



A 



nmt 



1, for t = m + n, \m — n\{y^ 0) 

2, for t = n + m — 2, ■ ■ ■ , \n — m\ + 2 
0, otherwise 



and the quantum numbers {Ij, J"} label the state beyond the ground state. The 
densities corresponding to charge rapidity and isospin rapidity on the real axis for which 
the "omitted k, A values" must be taken into account. By introducing the density of 
holes for charge p^{k) and isospin n-string o"^(A), we have 

1 1 



p + p 



TT 



L 



Ku2{k)+ / K,{k\k')p{k')dk' 



-E/ Kn/2{k\X)a4X)dX 



Kn/2{X) + J Kn/2{\\k)p{k)dk 
Y^^mnt I Kt,^{X\X')a^{\')d\' 



(32) 



where the integration limits are [0, oo]. In terms of the distribution functions of charge 
and isospin rapidities, the energy per length has the form Ek/L = J k'^p{k)dk, the total 
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number of down isospins is M/L = J2n^ I '^n{X)dX and the particle density of the 
system is D = N/L = f p{k)dk. 

If we consider the energy arising from the external field Vt which is the Rabi field 
in two-component BEC experiments, the internal energy of the system is 

E/L= f{k^- n)p{k)dk + ^ 2rifi f a„c/A (33) 

For a given p{k), p^{k), o"„(A) and o"^(A), the entropy is of the form 
S/L = J[{p + p^) ln(p + p^) - p In p - p^ In p'^jdk 



+ 5Z / 1^"^" + "^^^ ^^(^'^ + ^ri) - In (T„ - In a^jrfA 



(34) 



where the Boltzmann constant is set to unity. 

At finite temperature, we should minimize the free energy F = E — TS — pN 
where p is the chemical potential and 5* is the entropy of the system. Making use of the 
relations derived from eqs. (|32|) 

Spf'= -Sp+ I Ki{k\k')5pdk' -J2 I Kn/2{k\X)5and\ 

8ai = -Scr^+ I K^/2{X\k)Spdk-J2A^nt I Kt/2{X\X')SaM')dX' (35) 

mt 

we obtain the following conditions from the minimum condition 6F = 0, 

e{k) =e-n-p-T J Ki{k\k') ln(l + e~<'''^/^)dk' 

-TJ2 I i^„/2(A:|A)ln[l + e-^"W/^]rfA 

C„(A) =2nn + T j Kn/2{k\X) ln[l + e~'^''^/^]dk 

+ Tj2Amnt I K„2{X\X') ln[l + e-^-(^')/^]dA' (36) 

mt 

where we have written 
p^{k) 



pik) 



cr„(A) 

When T — > 0, eqs. (BB) become 



£o(fc) = e-n-p + j Ki{k\k')£o{k')dk' 



{31 



where the "ferromagnetic" ground state is under consideration. Clearly, the integral 
equation gives the solution for the dressed energy from which the ground-state 
energy can be given in terms of Sq 

Eo/L = - £o{k)dk (39) 

vr In 
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The solution of eq. (pS]) defines the energy bands. And the Fermi surface is determined 

by 

Soikp) = (40) 

since the ground-state configuration corresponds to the case that all states of So{k) < 
are fully filled. The bare energy ^^q^"* is the zero order term of eq. ( |5^ ) 

4°^ = e-n-fi (41) 

which is valid in the strong coupling limit 

Equations (^) can be solved by iteration. The coupled equations ( ^2|) of density 
distribution of charge and isospin are then a Fredholm type 

p(l + e'/T) = - - yK./^ik) + [ KMk')p{k')dk' 
n L ' J 

K^,MX)arXX)d\ 



O'r, 



[1 + e^"/^) = ^K„/2(A) + j K^/2{\\k)p{k)dk 

- / ^t/2(A|A0a„(A')ciA' (42) 



Finally, we obtain the Helmholtz free energy F = E — TS 



F = fxN-T j ln[l + e-^/^] - Ki/2ik) 
ln[l + e-'^"/^]ir„/2(A)(iA 



dk 



(43) 



and the pressure 



P = -% = - I ln[l + e-'/^]dk, (44) 



dL TT 

which is formally the same as Yang and Yang's expression but the equations determining 
e and ( are different. Consequently, the partition function is given by Z = e~^/^. The 
thermodynamic functions, such as partition function Z, free energy F, are of importance 
for a thermal system. Given either of them, one is able to calculate all thermodynamic 
properties for the system in principle. However, the eqs. (^) are so complicated that it 
is impossible to obtain the explicit form of e and ( for general case. Moreover, we can 
also obtain some plausible results for some special cases in the next section. 



6. Special cases 

In general, the free energy of our model should be calculated using formula (^), where 
e{k) and Cn(A) are determined from eqs. (|36|) . Eqs. (36) can be solved by numerical 
method via iteration. In this section, we will only discuss the case of strong coupling at 
low temperature. 
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When c oo, Kn{k\k') = and Kn{k) = 0, then we have 
e = A;^ — — /i 
and the free energy becomes 



F/L = ^D-^ j ln[l + e" 



which can be solved by integration by part [p!q] , 



F/L = f,D--f\ 



(45) 
(46) 

(47) 



24/iV2, 

where the external field is set to zero. 

We can not derive the specific heat directly from the afore-obtained free energy 
because the chemical potential is a function of temperature. From eqs. (^2]), the density 
of charge rapidity has the form 

1 1 

^1 + (,(k'^-n-f,)/T 



P 



(4^ 



Clearly, at zero temperature, the quasi- Fermi surface is just the square root of the 
chemical potential, so we have no = tt^D^ which denotes the chemical potential at zero 
temperature. At low temperature, however, it is determined by 



we have 



D 



1 + e('='-'^)/^ 



dk 



1 + 



TT 



24/^2 



(49) 



(50) 



The second term in the bracket is small at low temperature, so we can replace /x by /xq- 
The equation becomes 



= Po 



1 



then the free energy becomes 



F/L = fioD 



1 



12^1 



2 3/2 

in 



1 



(51) 



(52) 



Since by thermodynamics 5* = —dF/dT and = TdS/dT, we find the specific 
heat at low temperature is Fermi-liquid like 



S 



T 



(53) 



It is the same as the result of the one-component case, since for the strong coupling limit 
the isospin and charge are decoupled, the contribution of isospinon to the free energy 
vanishes. In figj^, the finite-size energy gap of holon-isospinon in strong coupling limit 
tends to zero. 
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7. Conclusions 

In this paper, we have solved a system of one dimensional trapped SU(2) bosons with 5- 
function interaction by means of the Bethe-ansatz method. On the basis of Bethe-ansatz 
equations we first discussed the ground state of the system and found that the ground 
state is an isospin "ferromagnetic" state which differs greatly from the spin-1/2 fermion 
systems. We studied the low excitation states by both numerical and analytic methods. 
It was shown that there are three elementary excitation modes, and the holon-antiholon 
and holon-isospinon excitations are gapless for large systems. For finite system, we not 
only plotted some excitation spectra but also plotted the dependence of finite-size energy 
gap on the inter-particle interaction. The thermodynamics of the system were studied 
by using the thermodynamic Bethe ansatz |jl3[. For strong coupling we found that the 



system exhibits the Fermi-liquid behavior, i.e. the specific heat is a linear function of T 
at low temperature. 
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Figure 1. The density of state in |fc|-space for the ground state (left). The distribution 

changes from a histogram to a narrow peak gradually for the coupling from strong to 
weak. The density of state in |fc|-spacc in the presence of one isospin rapidity at 
Ji = 30 where a rift evident appeared (right). The figure is plotted ior N = L = 100 
and c= 10,1,0.1,0.01. 




Figure 2. The ground state energy versus the coupling for different densities 
D = 1.0,0.75,0.5. 




Figure 3. The holon-antiholon excitation spectrum calculated for = L = 40 and 
c = 10 (left) and c = 1 (right). 
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Figure 4. The holon-isospinon excitation spectrum calculated ioi N = L = 20, c = 10 
(left) and c= 1 (right). 




I ' ' 1 

-10 -5 5 10 

interaction: log[c] 



Figure 5. The finite-size energy gap versus interaction for N — L — 50, obviously, 
c = corresponds to 0.02=1/L 




Figure 6. Finite-size effects: Isospin gap as a function of 100/L for D = N/L = 1, 
plotted for various inter-particle interaction 



